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Abstract 

The phase separation and vortex states in two-component Bose-Einstein con- 
densate consisting of \F = 1,771/ = —1 > and |2, 1 > internal spin states of 
S7 Rb atoms are considered in the framework of Thomas- Fermi approximation. 
It is shown that in nonrotating system the atoms in the state |1, — 1 > form 
a shell about the atoms in the state \2, 1 >. The critical angular velocity for 
each state is calculated. These velocities depend drastically on the relative 
concentrations of the components, the critical angular velocity of the outer 
component being less than the angular velocity of the inner one. It is shown 
that the atoms in the |1, — 1 > state can form a rotating ring about the resting 
core of the atoms in the state |2, 1 >. 

PACS:03.75.Fi; 05.30.Jp; 32.80.Pj. 

The realization of Bose-Einstein Condensation (BEC) in dilute atomic gases offers new 
opportunities for studying quantum degenerate fluids j|J|. These condensates, which contain 
thousands of atoms confined to microscale clouds, have similarities to superfluidity and laser, 
and provide new testing ground for many body physics. 

The modern theoretical description of dilute BEC originates from Bogoliubov's semi- 
nal 1947 paper where he showed that weak repulsive interaction qualitatively change the 



excitation spectra from quadratic free particle form to a linear phonon-like structure. To 
describe the trapped condensates at T = one can use the Gross-Pitaevskii (GP) (nonlinear 
Schrodinger) equation for the condensate wave function 0. This equation appears as the 
generalization of the Bogoliubov theory for the inhomogeneous phase. It was widely used 
to discuss the ground state properties and collective excitations in BEC. 

Bulk superffuids are distinguished from normal fluids by their ability to support dissi- 
pationless flow. This ability is closely related to the existence of stable quantized vortices. 
Such vortices have been widely studied in superfluid 4 He. Recently clear experimental ev- 
idences of the existence of a vortex in trapped BEC |3|-[5| was reported. Unlike superfluid 
helium, the trapping potential makes alkali BEC nonuniform. Theoretical work has been 
concentrated on the critical angular velocity of the vortex creation, the collective excitations 
of BEC in the presence of the vortex and considerations of stability of the vortex |J- [JT3j . 



It was shown that unlike superfluid helium where the vortex is locally stable, for weakly in- 
teracting gases stable quantized vortices exist only in a driven system, and become unstable 
without imposed rotation. So this system cannot be considered a superfluid P JIQyiT . 



In the study of quantum fluids, the most interesting behavior has been found in the 
behavior of fluid mixtures. At the present time, two experimental groups have observed 



trapped multiple condensates as realized in a magnetic trap in rubidium |L4| and in an 



optical trap in sodium fL5 |. In these experiments, the spatial separation of condensates has 
been observed. One can distinguish two types of spatial separation: (a) potential separation, 
caused by external trapping potentials; (b) phase separation, which can occur in the absence 
of external potentials due to the interaction between two condensates. It is the latter type 
of phase separation that we consider in this paper. This type of phase separation has been 
observed in the experiments on the simultaneously trapped condensates, consisting of the 
87 Rb atoms in the |2, 1 > and 1 1 , — 1 > spin states (states 2 and 1 correspondingly) ||14|| . 
In this case the intraspecies and interspecies scattering lengths denoted correspondingly as 
an, a 2 2, cl\2 are in the proportion an : a 12 : a 22 = 1-03 : 1 : 0.97 with the average of the three 
being 55(3)A |l||L6 . 



In this paper we consider the behavior of binary mixture of Bose-Einstein condensates 
of alkali atoms. We calculate the critical angular velocity needed to create stable vorticies 
in either component in the rotating frame. This quantity is of crucial importance in view of 
the experimental possibility of creating vortices by rotation of the confining trap ||. 

The physics of interpenetrating Bose fluids is very rich and far from a complete under- 
standing. While properties of a rotating single component Bose fluids are much studied ||- 
l3fl , the rotating mixture of alkali atoms presents a new frontier, that is essentially virgin 



territory. 

In order to derive analytic results, some approximations must be used. A commonly used 
one is the Thomas-Fermi Approximation (TEA) , which ignores the kinetic energy terms. It 
has been shown that in the case of one component condensates the TFA results agree well 
with the numerical calculations for large particle numbers, except for a small region near 
the boundary of the condensate [7||§. In fact, even for a small number of particles the TFA 
still usually gives qualitatively correct results. The TFA provides an excellent starting point 
of study. However, the TFA should not be relied upon when a quantitative comparison 
of experiment and theory is important. In this case a numerical approach based on the 
Monte-Carlo simulation becomes necessary. 

Let us first consider the phase separation in binary mixture without rotation. 

In the case of a two-species condensate, letting ^(r) (i = 1,2) be the wave function 
of species i with particle number iVj, we can write the two coupled nonlinear Schrodinger 
(Gross-Pitaevskii) equations as: 

h 2 l 

-VVi(r) + -mi^ix 2 + y 2 + AV)^(r) - ^i(r) + 



2m l r ' 7 2 

+ Gu|Vi(r)| 2 Vi(r) + G 12 |V 2 (r)| 2 ^i(r) = 0; (1) 

h 2 l 

V 2 V> 2 (r) + -m^\{x 2 + y 2 + X 2 z 2 )ip 2 (r) - + 



2m 2 2 

+ G 22 |V 2 (r)| 2 V 2 (r) + G 12 |Vi(r)| 2 ^ 2 (r) = 0. (2) 

Equations ([I]) and (0) were obtained by minimization of the energy functional of the 
trapped bosons of masses mi and m 2 given by: 



EfaM = / d 6 r 



^-|V^(r)| 2 + \m^\{x 2 + y 2 + AV)|^(r)| 2 + 



+ ^lW 2 (r)| 2 + \m 2 ul{x 2 + y 2 + AV)|^ 2 (r)| 2 + 

The chemical potentials /ii and /i 2 are determined by the relations / d 3 r\ipi\ 2 = iVj. 
The trap potential is approximated by an effective three-dimensional harmonic-oscillator 
potential well, which is cylindrically symmetric about z axis, A being the ratio of angular 
frequencies in the axial direction u zi to that in the transverse direction A = u Z i/ui. The 
experimental value of A is A = y/8. The interaction strengths, Gn, G 22 , G\ 2 are determined 
by the s-wave scattering lengths for binary collisions of like and unlike bosons: Ga = 
Anh 2 an / mi] G\ 2 = 2Trh 2 ai 2 /m, where m -1 = m{ 1 + raj 1 . 

Let us consider now the phase separation due to interaction between the two condensates. 
In this case 

\ m ^\ = \m 2 u 2 2 . (4) 

Let us simplify the equations by using dimensionless variables. We define the length 
scale 



a±=(—) 1/2 , (5) 



and the dimensionless variables 



r = a±r', (6) 
E = HoJiE', (7) 
ifc(r) = V / iV4^(r'). (8) 

The wave function tp'^r') is normalized to 1. In terms of these variables the Gross-Pitaevskii 
energy functional takes the form: 

E' = \j d\' [jVx| VVil 2 + N^x' 2 + y' 2 + AV 2 )K| 2 + iV 2 /3 2 | V'^| 2 + N 2 (x' 2 + y' 2 + AV 2 M| 2 + 

(9) 
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Here (3 2 = mi/m 2 = Lo\ju}\ and Ui = 8TranNi/a±. In deriving Eq. @ we used Eq. (|). Eqs. 
([[]) and (|2]) can be rewritten as: 



- vVi + (^ 2 + y' 2 + aV 2 )^ - a*xV4 + Vi + 
+ 



47rai 2 A^ 2 mi 2 

# 2 I A = 0; 



a± m 

?2v7'2„/./ i /„/2 , „,/2 , \2„/2\„/./ ,.'.,/ 



4 7 ra 12 iV 1 m 1 2 
H h/>i V>9 = 



(10) 



(11) 



a± m 
where n[ = 

In the TFA, Eqs. ([|), ((TOj) and (|Tl~D can be further simplified by omitting the kinetic 
energy. The phase segregated condensates do not overlap, so we can neglect the last terms 
in Eqs. (|9|), (|10"D and ( pi]) , obtaining from ( |10|) and (|Tl~D, in separate regions that they do 
not overlap, simple algebraic equations: 

1 



K(r')| 2 
l^(r')! 2 



"i 



£ - (p' 2 + AV 2 ' 



1 



M 2 /3 2 



a4 - (p' 2 + A 2 ,' 2 



(12) 
(13) 



Here p' 2 = x' 2 + y' 2 . From Eqs. (0) and ([□]) one can see that the condensate density has 
the ellipsoidal form. 

In the case of phase separation, the energy of the system can be written in the form 



E = E 1 +E 2 



(14) 



where 



E 2 



huiNi 
hu}iN 2 



(15) 
(16) 



In order to obtain Eqs. (]I5|)-(P^), Eqs. (p~^)- (p~3|) have been used. 

To investigate the phase separation in the mixture we first suppose that the condensate 
1 atoms form an ellipsoidal shell about the condensate 2 atoms (we will call this configu- 
ration as the configuration "a"). To determine the position of the boundary between the 



condensates, we use the condition of thermodynamic equilibrium |19 |: the pressures exerted 
by both condensates must be equal: 



The pressure is given by [ pOj] : 



Pi = P 2 - 



(17) 



(18) 



The condensate 2 has the form of the ellipsoid with long semiaxis q: 



p' 2 + AV 2 = q 2 . 



From Eqs. ([T2|)- (p~3|) and (P~7|)- (pT9|) one has the equation for q: 



(19) 



(20) 



where k = yj (a ll m 2 ) / (a 2 2m l ) . 

Chemical potentials fi[ and p 2 can be obtained using the normalization conditions 
/ d 3 r'\ip[\ 2 = J d 3 r'\ifj 2 \ 2 — 1 and are given by: 



(l - fg' 3 + §g' 5 



2/5 : 



// 2 



(M)3/2g'3 



15 5 



(21) 



(22) 



where q = J fi[q' and 



15 An,' 



2/5 



STY 



(23) 



From equations (^Tj)-(^) one can determine the chemical potentials fi[ and p 2 and the 
semiaxis of the phase boundary ellipsoid q as functions of Ni and N 2 . The energy of the 
configuration " a" E a = E al + E a2 is given by: 



E al = -h^N, - 



105 



3 5 y 7 



15 (K) 3 / 2 



£ a2 = ^iiV 2 <j fi' 2 - y ^^5/2 ^2 2 y - 2 /4/4^r + M 2 ^- 



/5 



17 s 

T 



(24) 
(25) 



Let us now consider the opposite case when the condensate 2 atoms form an ellipsoidal 
shell about the condensate 1 atoms (configuration "b"). In this case Eqs. (p0|)-(|25|) can be 
rewritten in the form: 



t4 



q\ = k(h% - q\) 



//\5/2 



W) 5/2 

5 J3 i 3 _/5 ' 
- 2?1 + 5?1 



(26) 
(27) 



2 (/i?)V2 



(28) 



£ 6 = 


E b i + E b 2 


-Em = 




Eb2 = 





15 (/4Q 3 ^ 2 [ mil o,,",, "1i I ,,"2<?i 7 
T(^V Ml Y ^iM 2 y + M 2 y 

_8_ _ / _ 2 gf \ 
105 I 3 5 yi 7 / 



15 y 2 Q 7/2 
4 /3 2 (^) 5 /2 



(29) 
(30) 

/x^i is the long 



Here /// and // 2 ' are the chemical potentials in the configuration "b", q\ = 
semiaxis of the boundary ellipsoid, Ej, is the energy of the configuration "b". 

To estimate which configuration is stable, one has to compare E a and Eb- Let us first 
consider the limiting cases n 2 = N2/N1 <C 1, and n\ — N1/N2 <C 1. 

In the former case n2 1 the approximate solution of Eqs. (EUI)- (22) has the form: 



/ii=/i°(l + gg), 



K Ho 



2 5 ^3 



6 



9o 



/4 



1 (l + (« - l)g 2 + g 3 ) . 



(31) 
(32) 
(33) 



where 



{2n2\ 1/3 



V 5/c" J 



(34) 



From Eqs. (|26|)-(j2"8"|) one has: 
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^ = ^(l + 3^o), 



1 



/4' = //° ( 1 - 2(« - l)p + 3 (6« 2 + 4k - l)pj 



or 



(35) 
(36) 
(37) 



where 



Po 



( 2n 2 \ 
U5k 2 



1/2 



J 



Using Eqs. fl3l"D -fl38D, it may be easily shown that 



1 1 — K 

AE = E a -E b = -^iJVi/i? 7i 2 . 

Z ft 



(38) 



(39) 



From Eq. (|39f) one can see that for k > 1, AE < 0, so the configuration "a" is stable. 
Let us now consider the case n\ <C 1. Approximate solution for the configuration "a" is 
given by: 



, „ 5(3 -2k) 2 (224k 2 - 630k + 225) s 
q = 1 - X ^ -X - - -X 



[ = ^1 + 2 



6k 72k 2 
[k - 1) (k 2 - Ak - 6) 



-x 



K 



3 k 



-X 



& = ^ 1 + V 



(40) 
(41) 
(42) 



Here 



'2K 2 ni' 
15 



1/2 



(43) 



In the configuration "b" solution has the form; 



(k-2/5) 3 5 4 

3k 6 



Mi = 



k — (k — l)y + Ky — 



2| 3 2(«-l)(5«-2)_ 



15k 



i4 = l& B A(i + i/ a 



(44) 

(45) 
(46) 



where 



/2«m^ 1/3 



(47) 



The energy difference is: 

AE = -hu^^n^l - k). (48) 
From Eqs. (B^) and (§51) it is seen that the configuration "a" has lower energy if k = 



«n^2)/(a22^i) > 1- For mi = m 2 this is consistent with the qualitative assertion and 
experimental observation that it is energetically favorable for the atoms with the larger 
scattering length to form a lower-density shell about the atoms with the smaller scattering 



length |J|2T 



To evaluate AE in general case it is worth first to estimate the energy of the phase 
boundary which arises due to the gradient terms omitted in the TFA. The surface energy 
per unit area, the surface tension, is defined as a = E s /S, where E s is the surface energy, 
and S is the interface area, a may be written in the form fll"Tyi8] : 



a 



nui l a 12 A 



( Ml n 2 iV 1 iV 2 ) 1 /4|^||^| (A r 1 |^|2 + NM\ 2 )^. (49) 



2\/2a 2 L \v /a ii a 22 / 

Taking into account that the surface area of the ellipsoid with the semiaxis a±q has the 
form: 

S = 2<KaW ( 1 + .\ log \ ) , (50) 

one can estimate the contribution of the surface energy E s = crS to the total energy of each 
configuration. To be specific, we will use the parameters corresponding to the experiments 
on 87 Rb atoms. In this case m\ = m 2 , a± = 2.4 x 10~ 4 cm, N = Ni + N 2 = 0.5 x 10 6 atoms. 

In Fig. 1(a) we show the energies of configurations "a" and "b" (including the surface 
energy) E a /(%jj\N) (solid line) and Eb/(hu)iN) (dashed line) as functions of log 10 ( n 2)- One 
can see that E a is always lower than E^. Fig. 1 (b) represents the difference AE = (E a — 
Efj) / (huiN) . The behavior of AE for small and large values of n 2 is well described by Eqs. 
(|39|) and flUf). Fig. 1 (c) illustrates the behavior of the surface energy as a function of n 2 . It 
should be noted that the surface energy is much smaller than the interaction energy because 
the scattering lengths have very close values (see Eq. (|49|) ). 

9 



Let us now consider a trap rotating with frequency along the z-axis. 

For vortex excitation with angular momentum hlj , the condensate wave function is given 

by 

^.(r) = |^ 3 (r)|e^. (51) 
In a rotating frame the energy functional of the system is 

E rot (h,l 2 ) = E^iM + / tfrtyl+^iKld^+ifa). (52) 



After substituting the wave function for the vortex excitation (|5l| ) in Eq. (|52|), the 
effective confinement potential for the bosons becomes l 2 h 2 /2m lP 2 + l 2 h 2 /2m 2 p 2 + V 1 + V 2 , 
where V{ = rriiUJi(p 2 + X 2 z 2 )/2 and p 2 = x 2 + y 2 . So within the TEA the density of the vortex 
state, in separate regions that they do not overlap, has the form: 

K(r')| 2 = - U{k) ~ (P' 2 + A 2 ,' 2 ) - 4) ! (53) 
ui V p' 2 



W(r')| 2 = — MM - (P" + A 2 ) - —J- . (54) 



The important new qualitative feature of a vortex in the TEA is the appearance of a small 
hole of radius Q oc If/piili), but the remainder of the condensate density is essentially 
unchanged. The fractional change in the chemical potentials caused by the vortex (/4(^) — 
p'ij/p'i can be shown to be small |||9| , of the order of 1/iV 4 / 5 . In the calculation of physical 
quantities involving the condensate density it is sufficient to retain the no-vortex density 
and simply cut off any divergent radial integrals at the appropriate core sizes £ 2 = l\/p'i 
or £| = (3 2 l 2 / 74- Note that using the unperturbed density for calculation of the vortex 
properties corresponds to the hydrodynamic limit. 

In the case of the phase segregated condensate, one finds from Eqs. ( pT| - |5*2| ) and flT5]-16) 
that the energy change due to the presence of the vortices AE = E rot (li, l 2 ) — E rot (0, 0) has 
the form: 



AE = AE Nl + 



N 2 



10 



= -fiwjiVi J d 3 r' 
+ -nuj x N 2 J d 3 r' 



> |2 



2n/i 



I |2 



'Z 2 2 /? 2 



V2 



I |2 



2f2Z 2 1 , / |2 



#2 



(55) 



In the hydrodynamic limit ^ is given by Eqs. ( O ) and ([13]). 

Let us consider the stable configuration "a". In the hydrodynamic limit the location of 
the phase boundary is given by Eq. (|i~9f ). From fl55f) one has: 



AE 



Ni 



/\3/2 



AE N2 

1 
2 



%UJ\N2 



W 2 \ 

2(^)5/2 



In 



2/4 

Zl 



3 , 
~2 q 



1 - ^g' 2 ) In 



2/ii<7' 



/1 



/4 - 



In 



2 V //i// 2 g 



Z 2 /3 



A*2 
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2fiZ 2 



(56) 



(57) 



The critical angular velocities required to produce the vortex states in each condensate 
can be determined from the conditions AEn ± < 0, AE^ 2 < and have the form: 



n 



Ni 



5Ji(a*'i 



M3/2 



Wl 2(/i?)V2 



In 



Zi 



4 N 
3 



3 , 
2 y 



1 



9 



/2 N 



In 



2/iW 
Zi 



9 



/2 N 



15Z 2 (M) V V 



Wi " 4(^)5/2 



/ - 1 - / /2 
»2 ~ 



In 



2V^i/^2? 



Z 2 /3 



A* 2 - 
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Let us consider the asymptotic behavior of the critical angular velocities (|58|) 
for N 2 <C Ni and iVi C 7V 2 . Using the approximate solutions (0)-(iD, we obtain: 

n Nl 5/, /, 2//? 1 



(58) 
(59) 
and P|) 



wi 2/i? V Zi 3 y 

fiAr 2 _ 3Z 2 /3 2 



Zi 



In 



' 2Ao x 



Z^ 2 
A*? 



In 



' 2/xggo 



4 

1 - -K 

5 



7 3 

K 

12 4 



(60) 
(61) 



From ( |60|) and (^T|) one can see that if 7x2 — ^ the critical angular velocity of the external 
condensate Vt^ 1 tends to that of the pure condensate with the scattering length an (see Eq. 
(26) in Ref. ||). The critical angular velocity of the inner condensate tends to infinity 
as n 2 —>■ 0. However, this consideration can not be applied to rapidly rotating gases with Q 
comparable to U\ where the form of the condensate depends on Q ]TJ) . 

In the opposite limit n\ — N1/N2 <C 1 the critical angular velocities can be written as 
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h 




2^ 4/5 




5Z 2 /3 6/5 




2$ 



ln 2^ 4/5 , h 



h 3k//§/? 4 / 5 



2/i°/? 4 / 5 ' 

(3 — «J In h 3 — 2ac 

'i 



X, (62) 



+^^-1^ (63) 



In deriving ( ]62| - |6"3"| ) we used the approximate solutions ( f40| - [43"| ) . Note that when m — > 0, 
the critical angular velocity ( |63"D has the same form as the critical velocity for the pure 
condensate with scattering length 022, HoP^ 5 being the chemical potential. 

Fig. 2 shows the critical angular velocities for the external (fijvi) and the inner (f2jv 2 ) 
condensates as functions of 712 = N2/N1 for l\ = I2 = 1. 

Using Eqs. (|52|) and (|55|) , one can find the vortex configurations which correspond to 
the energy minimum for a given angular velocity Q/u\. In Table 1 we represent angular 
momenta of condensates which correspond to the minimum of the total energy E tot /Uui of 
the system for different values of the angular velocity and n 2 . In calculation we use the 
parameters for the 87 Rb. E tot is calculated as a sum of E (Eq. (p^[), the surface energy E s 
(Eqs. (H-50)), AE Nl (Eq. and AE N , (Eq. (0)). 

To summarize, we have shown that in the case of 87 Rb condensate the atoms in the state 
1 form a shell about the atoms in the state 2, the critical angular velocity for each state 
being drastically dependent on the mutual concentrations. The critical angular velocity of 
the outer component is less than the angular velocity of the inner one. When the ratio of 
the number of state 2 atoms to the number of state 1 atoms is small enough the critical 
angular velocity of the inner state becomes very large: in the framework of the hydrodynamic 
approximation it turns out to be larger than the oscillator frequency characterizing the 
confining potential and to decrease smoothly with increasing number of atoms in the state 
2. From the Table 1 one can see that atoms in the state 1 can form a rotating ring around 
the resting core of the atoms in the state 2. 

This work was supported in part by NATO Grant No.PST.CLG.976038. V.N.R and 
E.E.T acknowledge the financial support from the Russian Science Foundation through the 
Grant No.98-02-16805. STC was partly supported by NASA under grant No. CRG8-1427. 
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TABLES 



TABLE I. 

fi/wi n 2 h h E tot 

0.1 1.0 1 1.296863 x 10 7 

0.1 10.0 2 1.287497 x 10 7 

0.15 1.0 2 1.295151 x 10 7 

0.2 1.0 3 1 1.292344 x 10 7 
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FIGURES 

FIG. 1. (a) The total energies of configurations "a" and "b" E a /(Tiuj\N) (solid line) and 
Eb/(htoiN) (dashed line) as functions of log 10 (n2). (b) The difference AE = (E a — E^) / (HujiN) as 
a function of log 10 (n2). (c) The surface energies as functions of log 10 (ri2). Solid line corresponds 
to the surface energy E sa for the configuration " a" , dashed line - to the surface energy E s f, for the 
configuration "b". 

FIG. 2. Critical angular velocities VLn 1 /uj\ and D,n 2 /u>2 as functions of log 10 (n2) for the 
configuration "a". Solid line corresponds to the outer condensate 1, dashed line - to the inner 
condensate 2. 
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